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Abstract 

We study the finite size corrections for tlie magnetization and the 
internal energy of the 2d Ising model in a magnetic field by using 
transfer matrix techniques. We compare these corrections with the 
functional form recently proposed by Delfino and LeClair-Mussardo 
for the finite temperature behaviour of one-point functions in inte- 
grable 2d quantum field theories. We find a perfect agreement be- 
tween theoretical expectations and numerical results. Assuming the 
proposed functional form as an input in our analysis we obtain a rel- 
evant improvement in the precision of the continuum limit estimates 
of both quantities. 
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1 Introduction 



Despite the impressive progress of these last years, 2d quantum field theo- 
ries (QFT) still provide many interesting open problems. One of these is 
the finite temperature behaviour of correlators. Besides its experimental and 
theoretical relevance, this issue is particularly important when these theo- 
ries are studied by means of numerical simulations. In fact choosing a finite 
temperature setting corresponds, in the euclidean formulation of the theory, 
to compactify the (imaginary) time direction on a circle whose circumfer- 
ence R coincides with the inverse temperature. Thus the finite temperature 
corrections become, in this framework, finite size corrections in the cylinder 
geometry which, as it is well known, play a crucial role in the extraction of 
continuum limit expectation values from numerical simulations. 

While for a generic 2d QFT finding a functional form for the finite tem- 
perature correction seems a hopeless task, recently two interesting propos- 
als by Delfino |I| and LeClair and Mussardo (based on a previous work 
by LeClair and collaborators [Q) appeared in the literature to address this 
problem in the simpler case of integrable models. In particular, in 0], 
the authors studied the finite temperature behaviour of one-point correla- 
tion functions. General arguments suggest that, outside the critical point, 
any one-point function evaluated on a finite size lattice should approach 
its infinite volume limit ($;)ij=oo with an exponential decay of the type 

~ ($i)R=oo + Ce"7 + . . . (1) 

^ being the correlation length (i.e. the inverse of the lowest mass of the the- 
ory). However with these considerations nothing can be said on the constant 
C and the higher order terms in the above equation. 

The main achievement of [|l], |^ was to show that the constant C is indeed 
universal and to predict its value. At the same time they were able to give 
an explicit expression for higher order corrections. 

The two proposals and have different theoretical starting points 
and indeed the predictions for these higher order corrections turn out to be 
different in the two approaches (for a critical comparison and a discussion of 
these differences see 0]). However, if one looks at the first few orders in a 
low temperature expansion (those which do not involve multi-particle form 
factors) the results coincide. 
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The aim of this paper is to test these proposals in a particular integrable 
model: the 2d Ising model in a magnetic field. This model has recently been 
the subject of several theoretical and numerical studies In particular, 

in the numerical analyses (performed using transfer matrix techniques or 
montecarlo simulations 0, §]) a great attention was devoted to the analysis of 
finite size correction. However, despite this careful treatment these correction 
were the major source of uncertainty in the final results. Our goal in this 
paper is thus twofold: first we shall use the numerical results for finite size 
lattices to test the proposals of [H, @]. Second, we shall see if, assuming the 
functional form for the finite size corrections proposed in we can improve 
the precision of the continuum limit results for the one -point correlation 
functions. 

This paper is organized as follows. In sect. 2 we shall briefly discuss the 
model, the observables and the proposal of [I], ^. Sect. 3 will be devoted to 
a discussion of our numerical test and to a comparison with Finally 
sect. 4 will be devoted to some concluding remarks. 



2 Ising model in a magnetic field 

The Ising model in a magnetic field is defined by the partition function 

Z = e^(E(„,™)'^"'^-+'f^En'^") (2) 

where the field variable an takes the values {±1}; n = (uq, rii) labels the sites 
of a square lattice size Lq and Li in the two directions and lattice spacing 
cQ. (n, m) denotes nearest neighbour sites on the lattice. In the following 
we shall treat asymmetrically the two directions. We shall denote no as the 
compactified "time" coordinate and rii as the space one. The number of 
sites of the lattice will be denoted by = LqLi. The lattice extent in the 
transverse direction will be denoted as i? = Lq. This length will play a major 
role in the following. In fact our goal will be to describe the R dependence 
of the expectation values of the spin and energy operators of the theory. 

^ Since the lattice spacing will play no role in the following we shall set a = 1 in the 
rest of the paper. 
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In order to select only the magnetic perturbation, the coupling (3 must 
be fixed to its critical value 

/3 = /3c = ^ log (72 + 1) = 0.4406868... 

by defining hi = (3cH we end up with 

Z{hl) = ^ e'^'=I^(n,m) <^™'^'"+''iEn'^" _ 

In the continuum limit the model is described by the action: 

A = Ao + h J £xa{x) . (4) 

where o"(x) is the perturbing operator and the perturbing field is the magnetic 
field h. 

In this paper we want to use results obtained in the continuum theory 
to describe the finite size scaling behaviour of the lattice model. In general 
this would require a precise definition of the relations between lattice and 
continuum quantities like for instance that between hi which appears in eq. (|^) 
and the perturbing field h. However, as we shall see below, we shall only be 
interested in adimensional ratios in which all these normalizations cancel out. 
Thus we shall neglect this problem in the following. The interested reader 
can find a careful treatment of this issue in 

2.1 Lattice operators 

As a first step of our analysis let us define the lattice analogous of the spin 
and energy operators of the continuum theory^. The simplest choices for 
these lattice analogous are 

^Strictly speaking these lattice operators do not correspond exactly to the continuum 
ones but are instead linear combinations of all possible relevant and irrelevant operators 
of the continuum theory with compatible symmetry properties with respect to the Z2 
symmetry of the model (odd for the spin operator and even for the energy one) . However 
near the critical point this linear combination is dominated by the corresponding relevant 
operator (cr for ai and e for e;) and the only remaining freedom will be a conversion 
constant relating the continuum and lattice versions of the two operators. As mentioned 
above we shall neglect in the following these normalization constants. 
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Spin operator 



(5) 



i.e. the operator which associates to each site of the lattice the value 
of the spin in that site. 

• Energy operator 




where the sum runs over the four nearest neighbour sites y of x. e?, 
represents a constant "bulk" term which we shall discuss below 0. 

The index I indicates that these are the lattice discretizations of the contin- 
uous operators. We shall denote in the following the normalized sum over all 
the sites of these operators simply as 

(^i = ^Y.Mx) = ^ ^ ei{x) . (7) 

In the following we shall be interested in the expectation value of these 
lattice operators. More precisely we shall be interested in: 

• Magnetization 

The magnetization per site M{hi) defined as 

^^^'^ ^ ^(^^))l^=/3. = ^(E^^) • (8) 

^ i 

which imphes 

M{hi) = {a{) . (9) 

■^As a matter of fact, for technical reasons, in our Transfer Matrix calculations we 
evaluated only the "time-like" part of the action i.e. '^no ni^(no,ni)^(no+i,ni)- While, 
obviously, this choice makes no difference in the thermodynamic limit, it might have some 
effect for finite values of R. We shall further discuss this point in sect. 3.1 below. 
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Internal Energy 

The internal energy density E{hi) defined as 

^(M = ^(E ^ncr^) ■ (10) 

(n,m) 

It is important to stress that in this case one also has to take into 
account a bulk analytic contribution (as it happens also for the free 
energy itself) Eb{hi) which is an even function of hi. Let us define 
eb = Eb{0). The value of Eb{0) can be easily evaluated (for instance by 
using Kramers- Wannier duality) to be = Let us define E{hi) = 

E{hi) — eb, we have 



E{hl) = J^{i: ^nCym) - . (11) 

{n,m) * 

Hence we have 

m) = {ei) . (12) 
In the following we shall in particular be interested in the ratios 

(13) 



{^i)r being the mean value on a lattice with transverse extent R of the 
lattice operator Since in the ratio all the normalization constants cancel 
out, we can identify the two ratios with the analogous ones evaluated in the 
continuum theory, with the continuum operators. 

{'^i)r {'^)r 



Notice that to perform this identification it is mandatory to eliminate the 
bulk correction as we did in eq.(|rT]) 

2.2 Critical behaviour 

The critical behaviour of magnetization and internal energy can be easily 
obtained by means of standard renormalization group methods. One finds 

M oc Ihl'^/y^-^ (15) 
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E oc I I ('^-f 



(16) 



where yt are the are the magnetic and thermal RG-exponents respectively. 

From the exact solution of the Ising model at the critical point we know 
that Vh = Y 2/t = 1- Inserting these values in the above expressions we 
find 

Moc\h\T5 (17) 

Eoc\h\T5 . (18) 

From the conformal field theory (CFT) description of the Ising model 
at the critical point it is possible to construct also the h dependence of the 
following terms in the two scaling functions ([17|) and (0). A detailed analysis 
can be found in 0. We shall make use of this result in the following. 

2.3 S-matrix results 

In 1989 A. Zamolodchikov suggested that the scaling limit of the Ising 
Model in a magnetic field could be described by a a scattering theory which 
contains eight different species of self-conjugated particles Aa, a = 1, . . . , 8 
with masses 

1712 = 2mi cos 

1713 = 2mi cos 

7714 = 2m2 cos 
ms = 2m2 cos 
mg = 2m2 cos 
m-j = 4m2 cos 
mg = 4m2 cos 

where mi{h) is the lowest mass of the theory, which coincides with the inverse 
of the (exponential) correlation length. After Zamolodchikov's paper several 



TC 
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1.6180339887..) mi, 
(1.9890437907..) mi, 
(2.4048671724..) mi, 
(2.9562952015..) mi, (19) 
(3.2183404585..) mi. 



30 

!^ cos — = (3.8911568233..) mi , 
5 30 ^ ^ 

lcos— = (4.7833861168..) mi , 
5 15 
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other interesting results were obtained, ranging from the exphcit values of 



various critical amplitudes |10, 11, 12| to the values of the overlap amplitudes 



of correlators |TB| , . All these predictions have been tested with numerical 



simulations both in the Id Ising quantum chain in the dilute IRF 
(Interaction Round a Face) model (which is another realization of the 
scaling Ising model in a magnetic field) and directly in the 2d lattice Ising 



model I, I, |T3, PI and in all full agreement between S-matrix 



predictions and numerical results was found. 

2.4 Delfino and LeClair-Mussardo proposals 

One of the most interesting features of the approach discussed in [0] is that, 
by using the form-factor technology the author was able to give a very explicit 
and compact expression for the first few orders of the finite size corrections 
both for {a) and (e) for the Ising model in a magnetic field. This result will 
be of great importance for our analysis. 

According to a generic one point function evaluated on a cylinder 
of transverse size R approaches exponentially its asymptotic value ($)_r=oo 
with the following law. 

= 1 + - E Mm^R) + 0(e-2-^) , (20) 

where mi,m2 and m^, are the first three masses (the only ones below the 
lowest pair creation threshold) of the Zamolodchikov's solution discussed 
above and Ko denotes the zeroth order modified Bessel function. 

The major result of fll was to show that the Af constants are indeed 
universal and can be evaluated exactly in the framework of the S-matrix 
description of the model. Their values in the case of the two operators in 
which we are interested here are given in tab. [l|. As mentioned in the 
introduction these results turn out to agree at this order with those obtained 
by LeClair and Mussardo in 0. 

3 Numerical test 

In order to test the above prediction we performed a numerical study of the 
finite size corrections by using the results of the transfer matrix analysis dis- 
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Table 1: The universal amplitudes entering the expansion ( |2^ ) for the Ising 
model in a magnetic held. 



$ 


a 


e 


Af 


-8.0999744.. 


-17.893304.. 


Al 


-21.206008.. 


-24.946727.. 


Af 


-32.045891.. 


-53.679951.. 



cussed in to which we refer for further details on the algorithm and on 
the raw data. We only recall here some general informations on the experi- 
ment. We studied the Ising model for the 19 different values of the external 
magnetic field listed in tab.§. For each choice of h we studied lattices with 
values of R ranging from 10 to 21. For each values of h and each lattice 
size we evaluated the mean value of the magnetization and internal energy 
(see however the footnote below eq.(^). We show in tab.|^ a typical sample 
of our data. An important ingredient of eq. (^OD are the values of the first 
three masses of the theory. In principle these masses could be evaluated di- 
rectly from the S-matrix solution, using the suitable normalization constants 
to match with the lattice discretization and then imposing the standard h 
dependence dictated by the RG analysis. However it is important to stress 
that at the level of precision of our analysis, these three masses show rather 
large corrections to scaling for the values of h that we study. To avoid large 
systematic errors in the analysis of the finite size correction it is mandatory to 
take into account these corrections (they are discussed in great detail in 
An alternative route, which turns out to be much simpler, is to evaluate these 
masses, when possible, directly from the transfer matrix calculation. For the 
present analysis we used this second route. We report for completeness these 
values in tab.^ and refer to P] for a detailed discussion of this table. 



3.1 Analysis of the data 

We performed a two steps analysis 
ll First of all, for each value of h we fitted the values of < a > and < e > 



as a function of R, according to eg . (|20|) : The results are reported in 



tab.^ Here are some comments on the fits: 
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Table 2: The Grst three masses. 



h 


1/mi 


l/m2 


l/ma 


0.20 


0.59778522553(1) 


0.37795775263(1) 


0.310888(1) 


0.19 


0.61388448719(1) 


0.38765653507(1) 


0.318578(1) 


0.18 


0.63134670477(1) 


0.39818995529(1) 


0.326940(1) 


0.17 


0.65037325706(1) 


0.40968266918(1) 


0.336077(2) 


0.16 


0.67120940172(1) 


0.42228634593(5) 


0.346115(3) 


0.15 


0.69415734924(1) 


0.43618773124(1) 


0.357209(3) 


0.14 


0.71959442645(1) 


0.45161985381(4) 


0.369548(4) 


0.13 


0.74799884641(1) 


0.4688779288(2) 


0.38338(1) 


0.12 


0.77998715416(1) 


0.488342470(1) 


0.3990(1) 


0.11 


0.81637015277(1) 


0.510513817(1) 


0.4168(1) 


0.10 


0.85823913569(5) 


0.5360654(1) 


0.4374(5) 


0.09 


0.9071039295(1) 


0.5659287(6) 


0.4624(5) 


0.08 


0.965123997(1) 


0.60144(1) 


0.492(1) 


0.075 


0.998514180(1) 


0.62189(1) 


0.508(1) 


0.066103019026467 


1.067300500(2) 


0.66405(5) 


0.543(1) 


0.055085849188723 


1.17524158(3) 


0.7305(1) 




0.044068679350978 


1.322589(6) 


0.82(1) 




0.033051509513233 


1.54057(2) 






0.022034339675489 


1.91(1) 







• We always kept the asymptotic value of < a > and < e > as a 
free parameter. 

• The error of the input data was always of the order of 10~^^. We 
started by fitting all the values of R and then eliminated them 
one by one, starting form the smallest ones, until we reached a 
reduced equal or smaller than 1. Then we accepted the result 
of the fit and stopped the analysis. These acceptable x^'s could 
be achieved for h > 0.075 keeping only the first term in eg. (|20|) . 
For h = 0.066.., 0.055.., 0.044.. we had to introduce also the second 
mass, and finally for h = 0.033.., 0.022.. all the three terms were 
needed^. The values of the masses which could not be directly 

''^This trend is due to the fact that we studied, for each value of h, the same range of 
values of R while mi decreases as h goes to 0. This means that the argument of the first 
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Table 3: Magnetization and Internal Energy at h = 0.075 



R 


M 


E 


10 


0.81 


S4069072534 





122792065846 


11 


0.8f 


S4096308760 





122797848782 


12 


0.81 


B4105802401 





122799890318 


13 


0.81 


B4109130761 





122800613003 


14 


0.81 


B4110302840 





122800869443 


15 


0.81 


S4110717053 





122800960633 


16 


O.Sf 


S4110863864 





122800993122 


17 


0.8f 


S4110916027 





122801004716 


18 


0.81 


S4110934601 





122801008861 


19 


0.81 


B4110941226 





122801010344 


20 


0.81 


B4110943594 





122801010875 


21 


0.8J 


B4110944441 





122801023168 



evaluated via transfer matrix (the empty spaces in tab.0) were 
evaluated using the S-matrix results and keeping into account the 
proper corrections to scaling. 

• For all the values of h the asymptotic values that we obtained 
for < (7 > and < e > from these fits are compatible within the 
errors with those quoted in p. In general the present estimates 
are more precise. This gain in precision becomes larger and larger 
as h decreases. For instance for h = 0.075 we found 

E = 0.122801011162(5) M = 0.884110944919(2) (21) 

to be compared with the values reported in p 

E = 0.1228010112(1) M = 0.88411094491(1) (22) 

with an improvement of more than one order of magnitude in 
precision. 

Let us stress that both this improvement and the fact that all our 
results agree within the errors (which means typically an agree- 
ment within 10 significative digits as in eq.s(^Tj) and ( |2^ ) ) with 

Bessel function in eq.(|2C|) (i.e. miR) decreases as /i ^ thus allowing to detect, within 
the errors of our data also higher order terms of the equation. 
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those of ref are rather impressive evidences of the correctness 
of eq. (pO[) p|. In fact, as it can be seen by looking at tab.^, the finite 
size corrections which eq. (^D|) is able to describe are much larger 
(five order of magnitude!) than the uncertainties of the data: the 
corrections are of the order of 10~^ while the precision of the input 
data is IQ-^^^ 

• For all the values of h we obtained stable and reliable values for 
Af. These values are reported in tab.^. The same was not true for 
Af. This is not strange. In all the similar analyses performed in 
it was impossible to have reliable estimates for the amplitudes of 
the subleading exponents due to the large systematic deviations 
induced by the uncertainties in the leading corrections. 

Looking at the data of tab.^ one can see that there are rather large 
corrections to scaling. This is true in particular for A\ (for which 
it is possible that the rather large magnitude of the correction is a 
consequence of our asymmetric definition for the internal energy, see 
footnote below eq. (^). In any case it is clear that in order to test 
the predicted values for A1 and A\ it is mandatory to discuss these 
deviations. To this end we performed a second level of analysis 

2] The aim of this second level of analysis is to use the data reported in 
tab.^ to reach the scaling limit values for A^ and A\. This problem is 
exactly the same that we had to face in to extract the continuum 
limit values of the critical amplitudes and thus we use here the same 
method developed and discussed in sect. 6 of to which we refer the 
interested reader for further details. 

We used as fitting functions 

Alihi) = Al{l + 6i,,|/i,|Tt + hAhi\^') (23) 



^The results of H were obtained with an iterative method (see sect. 5.1 of ref.fql and 



sect. 3. 1.2 of 17| for a discussion) based on the general assumptions on the finite size 
behaviour of one point functions summarized in eq.(P. In particular the discussion of 
sect. 3. 1.2 of [|l^ shows that the agreement with the present results is not a case but is due 
to the fact that the iterative method actually mimics the exact behaviour of eq.(pO|). It 
would be nice to see if this iterative algorithm keeps its predicting power even when the 
theory is not integrable and no exact prediction can be deduced from a S-matrix analysis. 
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Al{hi)=Al{l + hJhi\^o+b,^,\hi\'^o) . (24) 

Differently from the cases studied in keeping only two terms in the 
scaling function is enough due to the large errors in the input data. 

The final results are 

8.090 < \A1\ < 8.125 (25) 

17.0 < \A{\ < 18.1 (26) 

which perfectly agree with the predictions of |1|, 0] reported in tab.|l| 

Then (similarly to what we did in p) we tried to extract the dominant 
correction to scaling by assuming as fixed inputs the values for Af 
reported in tab.|l] and performing again the fits. The results are: 

0.56 < < 0.66 (27) 

- 0.96 < < -0.82 (28) 

As it also happens for the energy itself (see [0) the leading correction 
to scaling in the energy sector is rather large and negative in sign. 



3.2 An attempt to obtain A2 

It is impossible to extract the A^ by directly fitting. Moreover we cannot 
improve the situation by using the known continuum limit values of < $ > 
and Af due to the large correction to scaling whose uncertainty is of the 
same order of magnitude of the terms that we would like to observe. The 
only possible way out is to construct a combination of the input data which 
exactly eliminates < $ > and Af. This can be easily done by combining the 
values of $ measured at three different values of the lattice size Ri, R2, R3. 
The combination is the following: 

^ TT [^R^) - $(i?2)] Ai(l, 3) - [$(i?i) - ^(i^s)] Ai(l, 2) 
2 ($)h=oo A2(1,2)AK1,3)-A2(1,3)A,(1,2) ^ ^ 

where 

Ai{i,j) = [KoimiR,-KoimiRj)] . (30) 

The result obtained using as input data those for h = 0.022034... for the 
magnetization are reported in tab.^ The quality of the result is rather good 



12 



Table 4: Results of the Gts. 



hi 






0.20 


9.186(1) 


9.14(1) 


0.19 


9.1185(15) 


9.625(15) 


0.18 


9.0365(15) 


9.935(15) 


0.17 


8.9735(15) 


10.24(1) 


0.16 


8.912(1) 


10.55(1) 


0.15 


8.843(2) 


10.97(2) 


0.14 


8.770(4) 


11.44(3) 


0.13 


8.730(2) 


11.665(10) 


0.12 


8.660(5) 


12.035(10) 


0.11 


8.612(5) 


12.34(4) 


0.10 


8.562(2) 


12.64(2) 


0.09 


8.512(2) 


12.95(2) 


0.08 


8.460(2) 


13.35(5) 


0.075 


8.435(4) 


13.52(5) 


0.066103019026467 


8.389(2) 


13.83(5) 


0.055085849188723 


8.337(2) 


14.27(5) 


0.044068679350978 


8.288(2) 


14.65(5) 


0.033051509513233 


8.237(4) 


15.15(10) 


0.022034339675489 


8.206(20) 


15.68(15) 



and is compatible with the prediction for this constant reported in tab.|l[ 
Notice however that the quality of the results becomes worse and worse as 
h increases (they are already almost completely unstable at h = 0.055085...) 
moreover no result can be obtained for e even at h = 0.022034.... This can 
be easily understood. It due to the fact that the absolute value of < e > is 
much smaller than that of < a > and as a consequence the relative errors 
(which are those which generate the observed instabilities) are larger. We 
have checked that in our range of values of h higher contributions in the 
expansion of eq. (pOD essentially give no effect. 

4 Concluding remarks 

Let us briefly summarize our main results. 
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Table 5: Tentative estimate for 





i?2 




2 


15 


16 


17 


-22.3477682808347 


15 


16 


18 


-22.3151718797824 


15 


16 


19 


-22.2853020505359 


15 


16 


20 


-22.2593381627093 


15 


17 


18 


-22.2646875236694 


15 


17 


19 


-22.2261047816289 


15 


17 


20 


-22.1918745162029 


15 


18 


19 


-22.1652771489821 


15 


18 


20 


-22.1211891807013 


15 


19 


20 


-22.0505359976403 


16 


17 


18 


-22.1917546952317 


16 


17 


19 


-22.1416354502675 


16 


17 


20 


-22.0965339005886 


16 


18 


19 


-22.0640813462388 


16 


18 


20 


-22.0064100955253 


16 


19 


20 


-21.9155756522890 


17 


18 


19 


-21.9520757204636 


17 


18 


20 


-21.8778373894290 


17 


19 


20 


-21.7630532199712 


18 


19 


20 


-21.5973234560565 



• Our results nicely agree with the functional form for the finite size 
behaviour of the one-point functions proposed by Delfino and LeClair- 
Mussardo. We are also able to give a good estimate of the first correc- 
tion to scaling terms, both for (a) and for (e) and a rough but reliable 
estimate of the A2 for the magnetization. 

• Unfortunately we are not able to discriminate between the two pro- 
posals. This would require precisions which are definitely outside the 
range of our transfer matrix methods. 

• The asymptotic values for the energy and the magnetization obtained 
implementing eq.(pOD are more precise that those that quoted in of 
more than one order of magnitude. As the critical point is approached 
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this improvement in precision increases. 

• The agreement between the present results and those of 0, besides 
being a further strong evidence of the correctness of the proposal of [|T| 
and 0, also points out the power of the algorithm proposed in Q to 
deal with the finite size correction. This is an important result since 
this algorithm is not constrained to integrable theories 



It would be interesting to extend this analysis to other integrable mod- 
els which could offer examples in which it could be simpler to discriminate 
between the two proposals of and [0] and to compare them with other 
existing approaches (say, for instance, ||TP|1). It would also be important to 



extend the analysis to models which could be more easily accessible for com- 
parison with experimental results (see for instance the results of EO, ETI] on 



the Haldane- Gapped spin chains). At the same time it would be interesting 
to extend the present analysis to two-point functions, for which very precise 
numerical results exist for the 2d Ising model in a magnetic field ^ . This 
would be a rather important test since it has been recently claimed [p2i p3 



that the extension to two-point functions of the proposals and 0] could 
fail in the case of interacting theories. 
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